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A DYNAMICAL APPROACH TO THE LARGE-TIME BEHAVIOR OF 
SOLUTIONS TO WEAKLY COUPLED SYSTEMS OF 
HAMILTON-JACOBI EQUATIONS 



O " HIROYOSHI MITAKE AND HUNG V. TRAN 



Abstract. We investigate the large-time behavior of the value functions of the optimal 
control problems on the n-dimensional torus which appear in the dynamic programming 
for the system whose states are governed by random changes. From the point of view of 
the study on partial differential equations, it is equivalent to consider viscosity solutions 
of quasi-monotone weakly coupled systems of Hamilton- Jacobi equations. The large- 
time behavior of viscosity solutions of this problem has been recently studied by the 
authors and Camilli, Ley, Loreti, and Nguyen for some special cases, independently, but 
the general cases remain widely open. We establish a convergence result to asymptotic 
i-^h 1 solutions as time goes to infinity under rather general assumptions by using dynamical 

■ properties of value functions. 



>■ 1. Introduction and Main Result 

In this paper we deal with optimal control problems, or calculus of variations, which 
\q ■ appear in the dynamic programming for the system whose states are governed by random 

QO ■ changes. More precisely, we consider the minimizing problem: 

O 



Minimize E,- 



o 



Lu( s )(l{s), 7(5)) ds + 9v{- t ){l{-t)) 



'1.1) 



over all controls 7 G AC ([— t, 0]) with 7(0) = x for any fixed (x, t) G T n x [0, 00), where 
the Lagrangians Li(x,p) : T n x K" K are derived from the Fenchel-Legendre transforms 
of given Hamiltonians Hi and we denote by AC ([— t, 0]) the set of absolutely continuous 
functions on [—t, 0] with values in T n . The functions are given real-valued continuous 
functions on T™ for % — 1,2. Here Ej denotes the expectation of a process with u(0) = i, 
where v is a {l,2}-valued process which is a continuous-time Markov chain on (— 00, 0] 
(notice that time is reversed) such that for s < 0, As > 0, 

P(i/(s - As) = j I u(s) = i) = qAs + o(As) as As ->■ for i ^ j, (1.2) 

where q are given positive constants and o : [0, 00) — > [0, 00) is a function satisfying 
o{r)/r — > as r — > 0. We call the minimizing costs of (1.1) the value functions of optimal 
control problems (1.1). 
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The purpose of this paper is to investigate the large-time behavior of the value func- 
tions. From the point of view of partial differential equations it is equivalent to study 
that of viscosity solutions of quasi-monotone weakly coupled systems of Hamilton-Jacobi 
equations 

( ( Ul ) t + D Ul ) + ci(ui - u 2 ) = in T n x (0, oo), 

(C) \ {u 2 ) t + H 2 (x, Du 2 ) + c 2 {u 2 - ui) = in T n x (0, oo), 
[ui(x, 0) = gi(x) on T", 

where the Hamiltonians Hi(x, p) : T n x M n — > ffi. are given continuous functions for % — 1,2, 
which are assumed throughout the paper to satisfy the followings. 

(Al) The functions Hi are uniformly coercive in the x- variable, i.e., 

lim inf{Hi(x,p) \ x G M n , \p\ > r} = oo. 

r—>oo 

(A2) The functions p i— > Hi(x,p) are strictly convex for any x G T n . 

Here are real-valued unknown functions on T" x [0, oo) and {ui) t = dui/dt, Dui = 
(dui/dxi, . . . , dui/dx ri ) for i = 1,2, respectively. We are only dealing with viscosity 
solutions of Hamilton-Jacobi equations here and thus the term "viscosity" will be omitted 
henceforth. 

The existence and uniqueness results for weakly coupled systems (C) of Hamilton- 
Jacobi equations have been established by [6, 9]. In recent years, there have been many 
studies on the properties of viscosity solutions of weakly coupled systems of Hamilton- 
Jacobi equations. See [3, 11, 12, 4] for instance. In particular, the studies on large-time 
behaviors were done for some special cases by the authors [11], and Camilli, Ley, Loreti 
and Nguyen [4], independently. However, the general cases remain widely open and the 
techniques developed in [11, 4] are not applicable for general cases. The coupling terms 
cause serious difficulties, which will be explained in details later. 

Let us first recall the heuristic derivation of the large-time asymptotics for (C) discussed 
by the authors [11] for readers' convenience. We use the same notations as in [11]. For 
simplicity, we assume that c\ — c 2 — 1 henceforth. Formal asymptotic expansions of the 
solutions ui, u 2 of (C) are considered to be of the forms 

u\ (x, t) = a 01 (x)t + an (x) + a 21 (x)^ 1 + 
u 2 (x, t) = a 02 (x)t + a l2 (x) + a 22 (x)t~ l + . . . . 

as t — > oo. Then (C) becomes 

aoi(z) - a 2 i(x)t~ 2 + . . . + Hi{x, Da 01 (x)t + Da u (x) + Da^x)^ 1 + 
+ (a 01 (x) - a 02 (x))t + (a n (x) - a 12 (x)) + (a 21 (x) - a 22 (x))r l + ... 

and 

a 02 (x) - a 22 (x)t" 2 + . . . + H 2 (x, Da 02 (x)t + Da 12 (x) + Da^x)^ 1 + . . .) 

+ (a 02 (x) - a 01 (x))t + (a 12 (x) - a n (x)) + (a 22 (x) - a 21 (x))t~ l + . . . = 0. (1.4) 
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Sum up (1.3) and (1.4) to yield 

H x (x, Da 01 t + Da u + 0(l/t)) + H 2 (x, Da 02 t + Da 12 + 0(l/t)) + 0(1) = 

as t — >■ oo. Hence we formally get Da m = Da 02 = by the coercivity of Hi and H 2 . We 
next let t — >■ oo in (1.3), (1.4) to achieve that a 01 (x) = a 02 (x) = a for some constant 
a Q £ R, and 

f Hi(x, Dan(x)) + an(x) - a 12 (x) = -a , 

I H 2 (x,Dai 2 (x)) + (ii 2 (x) - an(x) = -a , 

in T n . It is then natural to study the ergodic problem 

f Hi(x, Dvi(x)) + vi - v 2 = c in T n , 
(E) < 

[ H 2 (x, Dv 2 {x)) + v 2 — v i = c in T n . 

We here seek for a triplet (vi,v 2 ,c) £ C(T n ) 2 x K. such that (^1,^2) is a solution of (E). 
If (fi,t>2,c) is such a triplet, we call (i>i,i>2) a pazr of ergodic functions and c an ergodic 
constant. It was proved in [3, 11] that there exists a unique constant c such that the 
ergodic problem (E) has continuous solutions (i>i,i>2)- 

Hence, our goal in this paper is to prove the following large-time asymptotics for (C). 

Theorem 1.1 (Main Result). Assume that (Al), (A2) hold. For any {gi,g 2 ) £ C{T n f 
there exists a solution (v 1 ,v 2 ,c) £ C(Y n ) 2 x IR of (E) such that 

Ui(x, t) + ct — Vi(x) — > uniformly on F as t — > 00 (1.5) 

for i = 1,2. 

In the last decade, the large time behavior of solutions of single Hamilton-Jacobi equa- 
tions, 

u t + H{x,Du) = in T n x (0,oo), (1.6) 

where H is coercive, has received much attention and general convergence results for 
solutions have been established. The first general result was discovered by Namah and 
Roquejoffre in [13] under the following additional assumptions: p t— > H(x,p) is convex, 
and 

H(x,p) > H(x,0) for all (x,p) eMxW n and max#(x,0) = 0, (1.7) 

M 

where M. is a smooth compact n-dimensional manifold without boundary. Then Fathi 
used dynamical system approach from weak KAM theory in [7] to establish the same type 
of convergence result, which requires uniform convexity (and smoothness) assumptions on 
H(x,-), i.e., D pp H(x,p) > al for all (x,p) £ M. x W 1 and a > but does not require 
the specific structure (1.7) of Hamiltonians. Afterwards Roquejoffre [14], Davini and 
Siconolfi in [5], Ishii in [8] refined and generalized the approach of Fathi and they studied 
the asymptotic problem for Hamilton-Jacobi equations on A4 or the whole n-dimensional 
Euclidean space. Besides, Barles and Souganidis [1] also obtained this type of results, for 
possibly non-convex Hamiltonians, by using a PDE method in the context of viscosity 
solutions. 
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In the previous paper [11], the authors could establish Theorem 1.1 only in two main 
specific cases. In the first case, we generalized the approach in [13] and obtain convergence 
result under additional assumptions similar to (1.7) (see also [4]). The second case is a 
generalization of [1] under the strong assumption that H\ = H 2 = H, where H satisfies 
similar assumptions as in [1]. We could not obtain Theorem 1.1 in its full generality 
because of the appearance of the coupling terms U\ — u 2 and u 2 — u\. 

In this paper we develop a dynamical approach to weakly coupled systems of Hamilton- 
Jacobi equations which is inspired by the works by Davini, Siconolfi [5] and Ishii [8], and 
establish Theorem 1.1 in its full generality. The results in [7, 14, 5] can be viewed as a 
particular case of Theorem 1.1 when H\ = H 2 , and g\ = g 2 . As we consider system (C), we 
need to take random switchings among the two states in (1.1) into account, which does 
never appear in the context of single Hamilton- Jacobi equations. The key ingredients 
in this approach consist of obtaining existence and stability results of extremal curves 
of (1.1). It is fairly straightforward to prove the existence of extremal curves by using 
techniques from calculus of variations. However, representation formulas (1.1) are implicit 
in some sense and prevent us from deriving a stability result (see Theorem 4.1). In order 
to over come this difficulty, we give more deterministic formulas for the value functions 
of (1.1) by explicit calculations in Theorem 2.4. By using the new formulas, which are 
more intuitive, we are able to derive Theorem 4.1, and hence large time behavior results. 

Let us call attention to the forthcoming paper [2] by Cagnetti, Gomes and the authors, 
which provides a completely new approach to the study of large time behaviors of both 
single and weakly coupled systems of Hamilton-Jacobi equations. A new and different 
proof of Theorem 1.1 is derived as well. 

The paper is organized as follows. In Section 2 we establish new representation formulas, 
which are more explicit and useful for our study here. We then derive the existence of 
extremal curves in Section 3, which is pretty standard in the theory of optimal control 
and calculus of variations. Section 4 concerns the study of stability of extremal curves. 
This section plays the key roles in this paper and allows us to overcome the technical 
difficulties coming from the coupling terms. See Remarks 4.4 and 4.5 for details. Section 
5 is devoted to the proof of Theorem 1.1. Finally, some lemmata concerning verifications 
of optimal control formulas for (C) in Section 2 are recorded in Appendix for readers' 
convenience. 

Acknowledgement. The authors are grateful to Professor Toshio Mikami for discussions 
which has been of help for them to come to Theorem 2.4. 



2. Preliminaries 

In this section, we establish new representation formulas, which give us a clearer intu- 
ition about the switching states of the systems. The new formulas allow us to perform 
deep studies on the extremal curves in Sections 3, and 4. 
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Lemma 2.1. Let v be a Markov process defined by (1.2) with c\ — C2 — 1 and set 
Pj(t) := P(i/(t) = j) for j G {1,2}. Then we have 

Pj(t) = 1/2 + e 2t (pj(0) - 1/2) for all t < 0. 

In particular, Pj(t) — > 1/2 as t ^ — oo /or an?/ j G {1, 2}. 

Proof Set ^(ax) := — (o"i(c?x) + ^((ix)), where Oj are the Dirac measures at the points 
j, for j G {1,2}, respectively. For any measurable function / and i = 1,2, by the Ito 
formula we have 

Ei[/(i/(f))] = /(t) +J o E i[J UXP) - /(*(*))) A*(#)] ^ 

= /(«') + = !) • (/( 2 ) - /C 1 )) + P W S ) = 2 ) • (/CO - A 2 ))} ^- 

If we take /(x) = l{j}(x) for j G {1,2}, where l^y is the characteristic function of {j}, 
then we have 

Pj (t)= Pj (0)+ [ -{p J ( s ).(-l) + (l- Pj (s))-l}ds=p 3 (0)+ [ (2 Pj (s)-l)ds. 
Jo Jo 

Therefore, we get Pj {t) = 1/2 + e 2t {pj{0) - 1/2) for all t < 0. □ 

A straightforward result of Lemma 2.1 is 
Lemma 2.2. Let 0j be any functions in C(T n ) for i = 1,2. We have 

HMt)(?)] = li 1 + z 2t )Hx) + ^(1 - e 2 %(x) 
for all x G T n , t < 0, and i = 1,2, where we take j so that {i,j} = {1,2}. 
Remark 2.3. In general if ci, c<i > are arbitrary constants, then we have 

UMt)( x )\ = —^—( c 3 + c^+^Ux) + -^—(1 - e^ +c ^)<P 3 (x) 

C\ + C2 C\ + C2 

for all x G T re , t < 0, and z = 1, 2, where we take j so that {2, j} = {1, 2}. 

It turns out that the value function of optimal control problems (1.1) can be written 
in more explicit forms without using continuous Markov chains as follows by using the 
Fubini theorem. 

Theorem 2.4. LetUi be the value functions defined by (1.1). Then we can write them as 

Ui (x, t) = inf { £ i(l + e 2 ')L i ( 7 (s), 7OO) da + i(l + e^^H)) 

+ /j(l-e 2s )L J ( 7 (4 7 ( S )) ( | S + |(l-^( 7 (-t)) I 7 GAC(H,0]), 7 (0) = *}. 

(2.1) 

Moreover, are uniformly continuous on T n x [0, 00) and t/ie £>azr (n 1; u 2 ) is the unique 
viscosity solution of (C) . 
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We call (1/2)(1 + e 2s ) and (1/2)(1 — e 2s ) for s < the weights corresponding to (C), 
which comes from the random switchings among the two states in (1.1). 

Proof. By Fubini's theorem and Lemma 2.2 we have 



E; 



L„( fl )(7(s),7(s)) ds + (7(-t)) 



E; 



^«/(a)(7( s )»7( s )) ds + E i[^(-*)(T(-*))] 
(1 + e 2s )L 4 ( 7 ( S ), 7(a)) d s + 1(1 + e - 2 ^( 7 (-t)) 



+ /° 5(1 - e 2s )LAl(s), i(s)) ds + 1(1 - e" 2 ')^( 7 (-t)) 

for any 7 G AC ([— t, 0]), which implies the equality (2.1). 

In Appendix we prove that are uniformly continuous on T n x [0, 00) and the pair 
M2) gives a solution of (C). In the previous paper [12], we showed that the pair (u\, u 2 ) 
defined by (1.1) solves (C) already. But we present it in a different way by using the new 
formula (2.1) itself to make the paper self-contained. □ 

Let (vi,V2,0) be a solution of (E). Without loss of generality, we may assume that the 
ergodic constant c = henceforth. We notice that Vi satisfies 



Vi(x) = inf \^i[J L v ( a )(7(s),7(s))ds + (7(-t))J | 

7 G AC ((-oo,0]) with 7(0) = x\, (2.2) 

where v is a {1,2}- valued process which is a continuous-time Markov chain satisfying 
(1.2) such that i/(0) = %. 

Proposition 2.5. Let (i>i, i>2, 0) be a subsolution o/(E). Then, 



vAx) < E 7 - 



^w(7(«)>7(s)) ds + (7(-t)) 



/or a// 1 > 0, 7 G AC ([-*, 0]) wift 7(0) = x. 

Lemma 2.6. Let t > 0, v t G W 1 '°°(T n ) /or i = 1,2 and j E AC ([— t, 0], T n ) iwtfi 
7(0) = x. VFe Ziaue u$ o 7 g AC ([— t, 0]) and there exists a function pi G L°°((—t, 0), R n ) 



Wi(x)=E 4 / p u (s)(s) -j{s) + ^2(v u (s) -Vj){j(s))ds + V v{ - t ){j(-t)) 

l J-t j=1 

Pi(s) G <9 c rj(7(s)) 



(2.3) 
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for a.e. s G (— t, 0) and i G {1,2}. Here d c Vi denotes the Clarke differential of Vi which 
is defined as 

d c Vi{x) = O co {Dvi(y) \ y G B(x,r), vi is differentiable at y} for x G T n , 

r>0 

where B(x,r) := {y G M n | \x — y\ < r}, and for A C IR n , coA denotes the closed convex 
hull of A. 

Proof. Fix any i G {1,2}. Let p G L7°°(IR n ) be a standard mollification kernel, i.e., p > 0, 
suppp C 5(0, 1) and f Rn p(x) dx = 1. Set p k (x) := k n p{kx) for fceN, and 

^ fc (z,t) = (*) := (p* * «0(7W) and p fc (z,t) = := 5(p fe * «i)(7(*)) 
for all t G (0, T). By the Ito formula for a jump process we have 



= E 5 



^(u(0),0)-^ k (u(-t),-t) 

r 2 



7 p k (u(s),s)^(s)ds+ [ J2(ip k (j,s)-i; k (v(s),s))ds 

J-t J-t j=1 

Note that ^ — >■ i>j(7(')) uniformly on [0,t] as — > oo and moreover passing to a 
subsequence if necessary, we may assume that for some Pi G L°°((0, T), M n ), p^ — ^ pj 
weakly star in i, 0)) as k — > oo, which implies (2.3). 

It remains to show that Pi(s) G 9 c Uj(7(s)) for a.e. s G (— £, 0). Since {pf}feeN is 
weakly convergent to pi in L 2 ((— t, 0), M n ), by the Mazur theorem, there is a sequence 
{??}fceN C L°°((-t,0),M n ) such that 

gf ->■ ft strongly in L 2 ((-t, 0), M n ) as fc ->■ oo, gf G co {p] | j > A;} (2.4) 

for all j G N. We may thus assume by its subsequence if necessary that 

Qi( s ) ~~ >" Pi( s ) f° r a - e - s e ( — ^0) as k oo. 

Now, noting that £>(p fc * t?)(x) = f yeB{Xjl/k) Pk{x - y)Dvi{y)dy for any x G T n and 
fc G N, we find that 

p k (s) G co{Dvi(y) | y G 5(7(5), l/k), Vi is differentiable at y} 
for any s G (— t, 0). Therefore, 

gf(s) G co^Pu^y) | y G 5(7(5), 1/&), is differentiable at y} 
for any s G (— t, 0). Since gf (s) — >■ Pi(s) for a.e. s G (— £, 0) as k — >■ 00, we get 

Pt(s) G O co~{Dvj(y) I y G 5(7(5), r), ^ is differentiable at y} = d c Vi(j(s)) 

r>0 

for a.e. s G (-t, 0). □ 
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Proof of Proposition 2.5. Let 7 G AC ([— t, 0]) with 7(0) = x and pi be the functions given 
by Lemma 2.6. In view of Lemma 2.6 we have 

f° 2 

E i / Pu(s)(s) -7(s) + ^2(v u (s) ~ Uj){j{s)) ds + V v{ - t )(^{-t)) 



VAX 



< E, 



< E; 



Hu( s )(l,Pu(s)) +•£'!/(«) (7> 7) -u j)(7)<& + «i/(-t)(7(-*)) (2-5) 



^(-)(7»T) ds + ^(-t)(7(-*)) 



3. Existence of Extremal Curves 



□ 



Let (fi,w 2 ,0) be a solution of (E). For any interval [a, b] C (— 00, 0], we denote by 
£(([a, 6], x, i, (vi, v 2 )) the set of all curves 7 6 AC ([a, 6]), which will be called an extremal 
curve on [a, b] such that 7(6) = x and for any [c, d] C [a, b], 



^i[vv(d){l{d))] =E i[y ^)(t( s )»7( s )) rfs + ^(c)(t(c)) 

with a continuous-time Markov chain v such that z/(0) = i and satisfies (1.2). 

Theorem 3.1. Let (1)1,1)2) be a solution of (E) . T/ien £((— 00, 0], a;, i, (v\, v 2 )) 7^ 0. 

In order to avoid technical difficulties we make the following additional assumptions in 
this section which are not necessary to get Theorem 3.1 and Theorem 1.1. We refer the 
readers to [8, Section 6] for the detail of general settings. 

(A3) Hi G C 2 (T n x W 1 ) and there exists 9 > such that D 2 pp Hi > 61 for i = 1, 2, where 

/ is the unit matrix of size n. 
(A4) There exists a constant C > such that 

^\p\ 2 ~C< Hi(x,p) < ^(\p\ 2 + 1) for x G T n , p G R n , % = 1, 2. 

Note that in this case we can easily see that Li G C 2 (T n x IR n ) are uniformly convex and 
satisfy 

7^\p\ 2 ~C < Li(x,p) <^(\p\ 2 + 1) for xeT n , peR n , z=l,2. (3.1) 
Lemma 3.2. Let (v\,v 2 ,0) be a solution of (E) . Then E([— 1,0], x,i,(vx,v 2 )) ^ $. 

Proof. By (2.2) there exists a sequence of curves {7^} C AC ([—1, 0]) with 7^(0) = x such 
that 

1 



Vi(x) + - > Ei 



^(«)(7fc(s)»7fc(s)) + ^(-i)(7fc(-l)) 



Since are bounded, we have 

-0 

E,: / L„( a )(7A.(s),7fc(s)) ds 



< C for some C > 0. 



(3.2) 
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Combining (3.2) and (3.1), we deduce that ||7/c||l 2 (-i,o) < M for some M > 0. For any 
— l<a<6<0, we have 



lk{b) - lk(a)\ < / \ik(s)\ds< / \j k (s)\ 2 ds 



1/2 



1 ds 



1/2 



< M\b-a 



1/2 



By the Arzela-Ascoli theorem and the weak compactness, by passing to a subsequence if 
necessary, {7*.} converges to 7 G AC ([—1,0]) uniformly, and {j k } converges weakly to 7 
in L 2 (-1,0). 

Now we prove that 



E, 



< lim inf Ej 

k— >oo 



^Ks)(7fe( s ))7fc( s )) ds ■ (3-3) 



This is a standard part in the theory of calculus of variations but let us present it here 
for the sake of clarity. The convexity of Lj gives us that 

Li(lk(s),ik{s)) > Li(~/ k (s), 7(s)) + D q Li(j k (s), j(s)) ■ (%(s) - 7(3)) 
= Li(-y k (s),j(s)) + [D q Li(y k (s), 7(3)) - D q Li{-y{s), j(s))] ■ (%(s) - 7(3)) 
+ 1)^(7(5), 7 (s)).(7 fc (s)- 7(5)). 

Since 7^ converges uniformly to 7, we employ the Lebesgue dominated convergence the- 
orem to get that 



lim Ej 

fc— >-oo 



Ly(a)(7fc(s),7(5)) (is 



Ei 



L„( s) (y(s),j(s))ds 



(3.4) 



We use (3.1) again to yield that 

\D q Li(x,q)\ < C(\q\ + 1) for x G T n , q G R n , z = 1,2. 

It it then straightforward by using the above and the Lebesgue dominated convergence 
theorem to see that 



lim Ej 

k— ¥00 



(D q L u(s) (y k (s),j(s)) - D q L u{s) (y(s),j(s))) ■ (<y k (s) -j(s))ds = (3.5) 



Besides, the weak convergence of {7*;} to 7 in L 2 (— 1,0) implies 



lim E, 

k— ¥00 



D q L u{s) (y(s),j(s)) ■ (%(s) -j(s))ds 



We combine (3.4), (3.5), and the above to get (3.3). Thus, 7 satisfies 



Vi(x) > E; 



L„(s){l(s),j{s)) ds + ^ ( _i)(7(-l)) 



(3.6) 
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On the other hand, for any — l<a<6<0, 

Vi(x) < E,; 



L v(s)(l(s),i(s)) ds + v u{b) (-f(b)) 
^iK{b)(j(b))} < Ej / L v{s) (j(s), j(s)) ds + v v{a) (7(a)) 

J a 

EiK( a )(7(a))] < Ej / L 1/ ( s )(7(s),7(s))ds + u v (-i)(7(-l)) 



The above inequalities together with (3.6) yield the conclusion that 7 G £([— 1, 0], x, i, (vi, V2)). 

□ 

Proof of Theorem 3.1. Fix x G T n and i G {1,2}. We define the sequence {7 fc }fc g N C 
AC ([—A;, —k + 1]) recursively as 7 fc G S([—k, —k + l],Xk-i, i,(vi, Vz)), where Xk := 7 fc (— k) 
and x = x. Define the curve 7 G AC ((—00, 0]) by 7(5) = 7 fc (s) for s G [—k, —k + 1] for 
k G N. Then it is clear to see that 7 G £((— 00, 0], x, i, (t>i,t>2))- □ 



4. Stability on the Extremal Curves 

In this section, we establish the following stability result, which plays a key role in the 
proof of Theorem 1.1. 

Theorem 4.1 (Scaling Result). Let (vi,v 2 ,0) be a solution of (E). For any r, T G 
(0,oo) with t < T such that t/(T — t) < 5q, where 5q appears in Lemma 4.3, and 
7 G £((—00, 0], x, i, (v\, V2)), we have 

Ui(x,T) -E i [« I/ (_ T )(7(-T),r)] 

<v t (x) - E^ ( _ T) ( 7 (-T))] + (1 + ^-^{JL-) (4.1) 

/or a fuction u : [0, 00) — > [0, 00) which is continuous and oo(0) = 0. 

Lemma 4.2. For any T > and 7 G £ ((— 00, 0], x, i, {v\, t> 2 )). There exists (^1,^2) £ 
L°°((-T,0),M n ) 2 suc/j inat 

Li( 7 (i), 7 (*)) + Hiijit),^)) = Pi (t) ■ j(t), 
Hi(r/(t),pi(t)) +Ui(7(*)) -«i(7(*)) = °' andPiW e cU( 7 (t)) 
/or a.e. t G (-T,0). 

Proo/. By Lemma 2.6 there exists (pi,p 2 ) G L°°((-T, 0), R™) 2 such that pi(i) G 9 c «t(7(*)) 
for a.e. t G (— T, 0) and satisfies (2.5) in the proof of Proposition 2.5. Also, note that by 
the convexity of Hi and the definition of Li, we have Hi(~f(t) , pi(t)) + Vi( r )(t)) — Vj{^f(t)) < 
and Hi(j(t),pi(t)) + L(^(t),^(t)) > Pi(t) ■ j(t) for a.e. t G (-T, 0) and z = 1,2. Since 
7 is an extremal curve, all inequalities above must become the equalities, which give the 
desired conclusion. □ 



LARGE TIME BEHAVIOR OF SOLUTIONS TO WEAKLY COUPLED SYSTEMS 



11 



Lemma 4.3. Let (wi,t> 2 ,0) be a solution o/(E). There exists So > such that for any 
e G [0, S ] and 7 G £((—00, 0], x, i, (vi, v 2 )) we have 



for a fuction u : [0, 00) — > [0, 00) which is continuous and w(0) = 0. 

Proof. Let (pi,p 2 ) be the pair of functions given by Lemma 4.2. We notice that 

iZi( 7 (t),Pi(t)) +H j (j(t), Pj (t)) = for a.e. t G (-00, 0] 

by Lemma 4.2. Set 

Q := {(x,j>i,j> 2 ) £T"x M 2n | ffxfopi) + # 2 (x,p 2 ) = 0}, 
S := {(x,gi,g 2 )| q% e D~Hi(x,pi) for some (x,p 1 ,p 2 ) G Q} 

and then Q and 5 are compact in T n x M 2n in view of the coercivity of iJj. We notice 
that (7(t),7(t),7(t)) G 5 for a.e. t G (-00, 0) and thus |-y(*) ] < M for some M > 0. We 
choose 5o G (0, 1) so that (x, (1 + e)j) G int (domLi fl domL 2 ) for all e G [0, 5q], where 
domLi := 6 P x R" h{x^) < 00}. 

By Lemma 4.2, 



Note that since Hi(x, •) are strictly convex, D q Li(x,^) exists, and is continuous on domLj. 

Due to the mean value theorem and (4.2), there exists 9 t G (0,1) and a fuction u> : 
[0, 00) — > [0, 00) which is continuous and u(0) = such that 



Remark 4.4. We notice that the result of Lemma 4.3 is different from the similar one for 
single equations (see [8, Lemma 7.2] for details). More precisely, the natural appearance 
of the coupling terms — e(vi — VjXtW) ma kes the analysis for weakly coupled systems 
more difficult. We could not proceed to establish large time behavior results in a crude 
way. It turns out that the weights (1/2) (1 + e 2t ) and (1/2)(1 — e 2t ) for t < are the key 
factors helping us overcome this difficulty as in the proof of Theorem 4.1 below. 



Proof of Theorem 4.1. Set e := r/(T-r) and T £ := T/(l+e). Notice that T = T £ + eT £ = 
T e + r. We have 

Ui (x, T) = ^(7(0), T) = Ui{y(0),T e + r) 



Li( 7 (t), (1 + e) 7 (t)) < (1 + e)L i ( 7 (t), 7 (t)) - e(v t - ^-)( 7 (t)) + eu{e) 



Li(j(t),j(t))= Pi (t)-j(t) - Hi(j{t), Pi {t)) 
D^it), y(t)) ■ j(t) + fa - VjWt)). 



(4.2) 




□ 



inf \ E; 
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Take 7 e £((— 00, 0], x, i, (fi,^)) and set r](s) := 7((1 + e)s) to derive that 



Mi (x,T)<E, / L i , (s) ( 7 ((l + £) S ),(l + £)7((l + £) S ))d S + w K _ Te) (7(-T),r) 



Make the change of variable t = (1 + e)s and use Lemma 2.2 to get 



Ui{x,T) <Ei 







T 1 + e 

x + e M/(l+e) 



^i/(t/(l+e))(7(*)> (! + £ )7<X>) d * + ^(-t £ )(7(-^), t) 

f° 1 _ g»/(l+«) 



(4.3) 



Li( 7 (t),(l + £) 7 (t))^ 



LMt),(l + e)j(t))dt 



_ T 2(1 + e) 
+ E t [^ ( _ Te) (7(-T),r)]. 

We use Lemma 4.3 in the above inequality to deduce 

Ui(x, T) - Ej [u I/ (_ Ts )(7(-T), r)] 

< y° ^(1 + e 2 '/^)I t ( 7 W,7W) + ^(1 - e 2 */(^))L,.( 7 (*), t(*)) * 



(4.4) 



+ 



1 + e 



,24/(14*)^ _ v .)( 7 (t)) dt + Tew(e). 



T 



We use the fact that u,- — is bounded in T n to derive that 



1 + e 



^/^\ Vj -Vi){r({t))dt 



T 



< Ce / e 2t/{1+e) dt < Ce. 



(4.5) 



(4.6) 



T 



Furthermore, for t < 0, |e 2 '/( 1+£ ) — e 2 '| < — 2tee 2 *^ 1+e \ This together with the facts that 
Mi are bounded and |-y(*) j < M imply Ej [u v ^- Te )(j(-T), r)] < Ej [m 1/( „ t) (7(-T), r)] 
and 

/ (e 2t/(1+£) - e 2 ')L fc ( 7 (t), 7 (t)) dt < -de [ te 2t/(1+e) dt < C 2 e, 

J-T J-T 

for k — 1, 2 and Ci, C 2 > independent of e. 
Summing up everything, we obtain 

Ui(x,T) -E i [tt I/ (_ T )(7(-T),r)] 



< 



rl x 



— T 



2 (1 + e 2i )L i ( 7 (t), lit)) + -(l-e^)L i ( T (t),7(t))j + Ce + Teo;(e) 



^(x)-E,[^ ( _ T) ( 7 (-T))] +C 



tT 1 t 

w = ), 



T-r T-r v T-r 



which is the desired conclusion. 



□ 



Remark 4.5. The new representation formula (2.1) with the weights (1/2) (1 + e 2t ) and 
(1/2) (1 — e 2t ) for t < appears naturally in both the statement and the proof of Theorem 
4.1 pointing out a major difference between single equations and weakly coupled systems. 
With new representation formula (2.1), we could explicitly calculate (4.4) and (4.5) and 
thus identify the main obstacle coming from the coupling term, the second last term in 
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(4.5). As mentioned in Remark 4.4, we could not estimate the coupling term in a crude 
way. For instance, in (4.3) we can easily see by Lemma 4.3 that 

"0 i 



< ^ 

< E, 



-L„(t/(i+s))(l(t), (1 + e)j{t)) dt 



T i + e 
o 



l + e 



~ «3-i/(t/(i+e))) dt + eTu(e) 



L, (t/(1+e)) ( 7 (t), 7(t)) dt + eTco(e) + CTe/(l + e) 



by using the fact that ||i>i — V2\\L°°(T n ) < C. But the last term in the above, which is of 
order 0{t) and does not vanish as s — > 0, is not enough to get the large-time asymptotics 
as we can see in the proof of Theorem 1.1. It turns out that the weights played an essential 
role here and helped us in establishing the key estimate (4.6) leading to the large time 
behavior result. 



5. The Proof of Convergence 
We define the functions u~. and (i = 1, 2) by 

u~.(x) := sup{<t>i(x)\ (4>i,(pj) is a subsolution of (E) with <pi < g i: <pj < gj}, 
u^x) := inf{-^j(x)| {jpi^j) is a solution of (E) with ^ > u~., ipj > 

where we take j so that {i,j} = {1,2}. Notice first that the pair (m 1; u 2 ) is a solution 
of (E) because of the convexity of the Hamiltonians Hi for i — 1,2. Let (T t )t>o be the 
semigroup defined by (T t gi)(x) := Ui(x,t), where (ui,n 2 ) is the unique solution of (C) 
with a given initial data (pi,^)- 

Lemma 5.1. We have (T t u~)(x) = mf s > t Ui(x,s), and u^x) = liminf^oo Ui(x, s) for 
x GT n audi = 1,2. 

Remark 5.2. In view of the uniform continuity of Ui on T n x [0, oo), we have liminf^oo Ui(x, 
liminf ^^u^x, t) := \im r ^oinf{ui(y, s) \ \x — y\ < r,s > 1/r} for all x G T n . 

Proof. Let Wi(x,t) := mi s > t Ui(x,s) for i — 1,2 and (x, t) G T n x [0,oo). Noting that 
{u~ u~) is a subsolution of (E) and u~ < gi, we get by the comparison principle 

u~(x) < (T s g~)(x) < (T s gi)(x) = Ui(x,s). 

We use the comparison principle again to deduce that 

(TfM~)(x) < (T t Ui)(x, s) = Ui(x,t + s) for all s > 0. 

Hence (T t u~.)(x) < mf s > t Ui(x, s) = Wi(x,t) for all (x,t) G T n x [0,oo). 

On the other hand, (wi,?/^) is a solution of (C) by the convexity of Hi. Moreover, Wi 
are increasing in the t- variable for % — 1, 2 by definition, which give us in addition that 
(wi(-,t),W2(-,t)) is a subsolution of (E) for all t > 0. In particular, (wi(-, 0), W2(-, 0)) is a 
subsolution of (E) with Wi < gi, which implies Wi < u~.. Thus, the comparison principle 
yields Wi{-,t) < (T t v,-)(-) on T n for all t > 0. 
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Next, note that 

liminf Ui(x, s) = lim (T t u~)(x) < lim (T t u^){x) = u^x) for i — 1, 2, x G TP 1 . 

Due to the fact that {T t u~_){x) = Wi(x,t) are increasing in the t-variable for % — 1,2, we 
see that (lim^oo u>i(-, t), lim^oo w 2 (-, t)) is a solution of (E). Therefore, by the definition 
of (m!,m 2 ) we deduce that lim^oo Wi(-, t) > Wj on T n . □ 

Proof of Theorem 1.1. Set Ui(x) := limsup^^ Ui(x, t). By Lemma (5.1) we have < U{ 
on TP 1 . We now prove that = Hi on T n for i = 1,2. 

Were the above statement false, there would exist a point igT" such that 

maxmaxfuj — uA (z) = (v,i — U;)(x) =: a > 0. 

Take 7 G £((— 00, 0], x, i, iux^u^)). We can choose a sequence {T m } C (0, 00) converging 
to 00 such that lim m _ ) . (X) Ui(x, T m ) = Ui(x) > u^x). Without loss of generality, we assume 
further that r y{— T m ) — > y G T n as m — > 00. 

Fix 5 > and choose r > large enough such that Ui(y, r) < u i (y)+a/2 and Uj(y, r) < 
Uj(y) + a/2. Notice that the first constraint requires r to be specific while the second 
constraint only requires r to be large enough. We apply Theorem 4.1 to get 

Ui {x,T m ) - {1(1 + e- 2T ™)u i { 1 {-T m ),T) + 1(1 - e- 2T -) Mj ( 7 (-T m ),r)} 
<n^) ~ {~(1 + e- 2T -h,( 7 (-T m )) + 1(1 - e- 2 ^( 7 (-X m ))} + (1 + 

I. Z Z J 1 m T 1 m 

for m large enough. Let m — > 00 in the above inequality to yield 

Ui(x) - -(Ui(y, r) + Uj (y, r)) < u^x) - -(u^y) + u^y)), 
which contradicts the choice of r. This finishes the proof. □ 

6. Appendix 

Let Ui be the value function associated with (1.1), or equivalently the function defined 
by the right hand side of (2.1). 

Proposition 6.1 (Dynamic Programming Principle). For any x G W l , < h < t and 

i = 1,2, we have 

Ui (x, t) = inf { J° 1(1 + e 2s )Z,( 7 (s), 7(5)) ds + 1(1 + e~ 2/l H (7 (-/*), t - h) 
+fj 1 {l-e 2s )L ] { 1 {s)^{s))ds+^ I T e AC ([-/*, 0]), 7(0) = x 



(6.i; 



Proof. We denote by tt>;(x, t, /i) the right hand side of (6.1). For any 7 G AC ([— t, 0]) with 
7(0) = x, set 7](s) = 7(5 + /i) for s G [— £ + h, 0]. Note that for s < 0, 

1(1 ± e 2 ^)) = 1(1 + e" 2 ") • 1(1 ± e 2 *) + 1(1 - e^) • 1(1 T e 2 *), 
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which actually comes from the memoryless property of Markov processes. We have 
i(l + e 2 -)L i ( 7 ( a ), 7W) + 5(1 + e- 2t )gM-t)) 
+ - e 2s )L,( 7 ( S ), 7 (s)) ds + 1(1 - e -«) ft .( 7 (-t)) 



o 



o 



(l + e^)L 4 ( 7 ( a ),7( a ))d«+ / -(1 - e^L^a), 7 (a)) da 



-(1 + e^LMs), V(s)) ds+-(l + e-^) 9l ( V (-t + h)) 

t+h * * 



+ 



t+h 



-(1 - e^foOO, d S + -(1 - e- a ^)^(»7(-* + h)) 



-(1 - e^M*), 77(a)) ds + -(1 - e -a(*-fc)) ft (,(_ t + h )) 



t+h 



+ 



t+h 



-(1 + e 2s )LMs), V(s)) ds + -(l + e-^) 9j (rj(-t + h)) 







1 



> / - (1 + OL^), 7(a)) ^ + g (1 + e- 2/l K(7(-/>), t " *) 
r ° 2 v- 3 2s )L i ( T ( S ), 7 (a))d S + ^ - 2ft 



+ / 1(1 - e 2s )L,( 7 (,), 7(a)) da + ±(1 - e"** )^( 7 (-^), t - h), 



which implies Ui(x, t) > Wi(x, t, h) for i = 1, 2. 

We also can prove the other inequalities by a similar way Thus, we omit the details. □ 

Proposition 6.2. The functions Ui are continuous on T n x [0, 00). 

Proof. We first prove that Ui are Lipschitz continuous on T n x [0, 00) under the additional 
assumption that 9i are Lipschitz continuous on T n . This additional requirement on 9i will 
be removed at the end of the proof. 

We may choose a constant M\ > so that Hi(x, Dgiix)) + (<& — gj)(x) < Mi for a.e. 
x G T n . It is clear that the function Vi(x, t) := Qi(x) — M x t on T n x [0, 00) is a subsolution 
of (C). 

By a similar argument to the proof of Proposition 2.5 we obtain 



Vi(x,t)<Mi / L iy(s) ( 7 (s),7(s))ds + ^(^)( 7 (-t),0) 



for all (x,t) G T n x [0, 00) and 7 G AC ([— £, 0]) with 7(0) = x, from which we get 
gi(x) - M\t < Ui{x,t) for all (x,t) G T n x [0, 00). 
It follows from (2.1) that 



Ui(x,t) <Ei / L v ( 8 )(x,0)ds + g v (- t )(x) < gi(x) + dt 
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for C\ : = max{Mi, maXi = i i 2;a;eT' 1 1-^(^0)1 + max xg T« \gi(x) — g2(x)\} , and all (x,t) G 
T n x [0, oo). Therefore we get 

\ui(x,t) -gi{x)\ < Cit for all (x,t) GT"x [0,oo). (6.2) 

Now, for any (x,t) G T™ x (0, oo) and h > 0, by Dynamic Programming Principle (6.1), 

r r° 1 1 

Ui (x,t + h) = M{ J -{l + e 2s )L l { 1 {s)^{s))ds + -{l + e- 2t )u l { 1 {-t),h) 

+ £i(l-e 2 %( 7 ( S ),7( S ))^ + ^(l-e-> J ( 7 H),/ i ) | 7 GAC(H,0]), 7 (0) = x). 

By (6.2), \ui{ r y{— t), h) — <?i(7(— t), h)\ < G\h for i — 1,2. Hence, we derive that 

i + /i) — t)| < C\h. (6.3) 

We next prove that Ui are Lipschitz continuous in x for i = 1,2. Fix x, y G TP with 
x y and £ > 0. In view of the coercivity of Hi, there exists a constant p > such that 
Li(x,£) < C for all x G T n and f G 5(0, p) (see [8, Proposition 2.1]). Set r := \x - y\/p 
and we first consider the case where r < t. Set rj(s) := y — sp(x — y)/\x — y\ for s G [— r, 0]. 
Note that t] G AC ([—t, 0]), 77(0) = ?/ and t)(—t) = x. By Dynamic Programming Principle 
(6.1), (6.2) and (6.3), 

Ui(y,t)<Ei J L v ( s )(r](s),r](s))ds + u u (- T )(r](-T),t-T) 

< Ct + -(1 + e~ 2T )ui(x, t-r) + -(l- e~ 2T ) Uj (x, t - r) 

< (C + 2Cit)r + Uj(x, t) < C|x -y\+ u^x, t), 

By symmetry we conclude \ui(x,t) — Ui(y,t)\ < C\x — y\, where C depends on t as 
calculated above. Notice that this is just a fairly crude estimate, but it is good enough 
for our presentation here. 

We consider the case where t < r. By (6.2), 

\ui(x,t) -Ui(y,t)\ < \ui(x,t) -gi(x)\ + |flfi(ar) - gi(y) \ + \gt(y) - Ui(y,t)\ 
< 2Ct + M\x-y\ < C\x-y\, 

where M := max i= i i2 ||^^j||L oo (T n )- Thus, we get \ui(x,t) — Ui(y,t)\ < C\x — y\ for all 
x,y G T n ,t > and z = 1,2. 

We finally remark that we can deduce the continuity of by using an approximation 
argument. We may choose a sequence {gf}keN of Lipschitz continuous functions so that 
~ 9i\\L°°{T n ) < 1 /k for all k G N. Let be functions defined by (2.1) with given 
g\. By comparison through the formulas for Ui and uf, we see that \m(x,t) — u^(x, t)\ < 
max T » \gi — gf\. Since u\ G C(T™ x [0, 00)) by the above argument for all k G N and 
converges uniformly to Uj on T n x [0, 00), we obtain Ui G C(T n x [0, 00)). □ 

Proof of Theorem 2.4. It is clear that {u\, «2)(-, 0) = (<7i,<72) on T n . We now prove that 
u\ is a subsolution of (C). Take a test function </> G C 1 (T n x (0, 00)) such that u\ — <fi has 
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a maximum at (xo,to) 6T"x (0, oo) and (u\ — <p)(xo,to) = 0. Take h > small enough. 
By Proposition 6.1, 



«i(x 0> t ) < 1° \{l + e 2s )L 1 ( 7 (s), 7(a)) + e^H^M), t - &) 

+ ^(1 - e 2s )L 2 ( 7 ( S ), 7(a)) * + ^(1 - e~ 2h )uM~h),t - h) 

for any 7 G AC([— /i, 0]) with 7(0) = Xq and 7 (0) = q G M n . We now use the fact 
u i(l(—h),t — h) < <f>(y(—h),to — h) to plug into the above to derive that 



0(7(O),to)-0(7HO,*o->O / 1 f" 1 , 



-X ./ ,2 (1 + C )Mt(s),7(s))^ 

1 /"° 1 1 - e~ 2h 

+ hj h 2 {1 ~ ^W^tOO) + ~^—( u 2 ~ ux)(y(-h),t - h). 

Sending h — > 0, we obtain 

<M^o,t ) + -D0(x o ,t o ) •<? < Lx(x ,q) + (w 2 -Ui)(x ,to) for all q G M", 

which implies (f> t (x ,t ) + Hx(x , D4>(x ,t )) + (u x - u 2 )(x ,t ) < 0. 

Next we prove that U\ is a supersolution of (C). Take a test function G C 1 (T n x (0, 00)) 
such that Mi — has a minimum at (xo, to) G T n x (0, 00) and (wi — 4>)(xo, to) = 0. Take 
h > small enough. By Proposition 6.1, there exists 7^ G AC ([— h, 0]) with 7/1(0) = xq 
such that 



/■o x 1 

Ui(z ,*o) + h 2 > -(1 + e 2s )L 1 ( 7 / i (s),7/ l (s)) + - (1 + e~ 2 >i ( 7A (-/i), to - &) 
J-/ t ^ ^ 

+ |° \{l ~ e 2s )L 2 ( lh (s),%(s)) ds + i(l - e^)n 2 ( 7 ,(-/i), to - fc)- 
We use ui( 7 ft(— /i), to — h) > 4>{lh{~ h), to — fo) in the above to yield 



0(7fc(O),to)-0(7fc(-/i),*o-/i) , ,^1 



+ />> x / -(l + e 2s )L 1 (7/ l ( S ),7^(a))^ 
-/ -(1- e 2s )L 2 ( T ,( S ), ^( a ))d a +_^(« 2 - til )( 7 (-/i) l t -/,). (6.4) 

J — h 
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On the other hand, 

<p(lh(0),t ) - <f){ lh {-h),t - h) 



1 



h 

o 

4>t(j h (s),to - s) + D4>(j h (s),t - s) ■ 7 ft (s) ds 



1 '° 



< r / Mlh(s),t - s)ds 



+ I J J (1 + ^{^OO, 7fc(*)) + #i(7fc00, ^(7h(s), *o - s))} da 

L 2 ( 7 ft( S ),7ft( S )) + F 2 (7ft( S ), J D0( 7 ft( S ),t o - s))} ds. (6.5) 



1 /- 1 



Combine (6.4), (6.5) and then send h — > to get 

(f>t(xo,t ) + ifi(x o ,-D0(x o ,t o )) + («i - M 2 )(x ,to) > 0. 
It is easy to see the uniform continuity of Ui due to the coercivity of Hamiltonians. □ 

References 

1. G. Barles and P. E. Souganidis, On the large time behavior of solutions of Hamilton- J 'acobi equations, 
SIAM J. Math. Anal. 31 (2000), no. 4, 925-939. 

2. F. Cagnetti, D. Gomes, H. Mitake and H. V. Tran, paper to appear. 

3. F. Cagnetti, D. Gomes and H. V. Tran, Adjoint methods for obstacle problems and weakly coupled 
systems of PDE, to appear in ESAIM Control Optim. Calc. Var. 

4. F. Camilli, O. Ley, P. Loreti and V. Nguyen, Large time behavior of weakly coupled systems of first- 
order Hamilton- J acobi equations, to appear in NoDEA. 

5. A. Davini and A. Siconolfi, A generalized dynamical approach to the large-time behavior of solutions 
of Hamilton- J acobi equations, SIAM J. Math. Anal. 38 (2006), no. 2, 478-502. 

6. H. Engler, S. M. Lenhart, Viscosity solutions for weakly coupled systems of Hamilton- Jacobi equations, 
Proc. London Math. Soc. (3) 63 (1991), no. 1, 212-240. 

7. A. Fathi, Sur la convergence du semi-groupe de Lax-Oleinik, C. R. Acad. Sci. Paris Ser. I Math. 327 
(1998), no. 3, 267-270. 

8. H. Ishii, Asymptotic solutions for large-time of Hamilton- Jacobi equations in Euclidean n space, Ann. 
Inst. H. Poincare Anal. Non Lineaire, 25 (2008), no 2, 231-266. 

9. H. Ishii and S. Koike, Viscosity solutions for monotone systems of second-order elliptic PDEs, Comm. 
Partial Differential Equations 16 (1991), no. 6-7, 1095-1128. 

10. P.-L., Lions, G. Papanicolaou and S. R. S. Varadhan, Homogenization of Hamilton-Jacobi equations, 
unpublished work (1987). 

11. H. Mitake, H. V. Tran, Remarks on the large-time behavior of viscosity solutions of quasi-monotone 
weakly coupled systems of Hamilton-Jacobi equations, Asymptot. Anal., 77 (2012), 43-70. 

12. H. Mitake, H. V. Tran, Homogenization of weakly coupled systems of Hamilton-Jacobi equations with 
fast switching rates, submitted. 

13. G. Namah and J.-M. Roqucjoffrc, Remarks on the long time behaviour of the solutions of Hamilton- 
Jacobi equations, Comm. Partial Differential Equations 24 (1999), no. 5-6, 883-893. 

14. J.-M. Roqucjoffrc, Convergence to steady states or periodic solutions in a class of Hamilton-Jacobi 
equations, J. Math. Pures Appl. (9) 80 (2001), no. 1, 85-104. 



LARGE TIME BEHAVIOR OF SOLUTIONS TO WEAKLY COUPLED SYSTEMS 



19 



(H. Mitake) Department of Applied Mathematics, Faculty of Science, Fukuoka Univer- 
sity, Fukuoka 814-0180, Japan 

E-mail address: mitake@math.sci.fukuoka-u.ac.jp 

(H. V. Tran) Department of Mathematics, The University of Chicago, 5734 S. University 
Avenue Chicago, Illinois 60637, USA 

E-mail address: tvhung@math.berkeley.edu 



